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PREDICTION OF ATTENUATION OF THE 28 GHz COMSTAR BEACON SIGNAL 
USING RADAR AND MEASURED RAIN DROP SPECTRA 


Abstract 

A program to measure the rain attenuation of the CC^TAR beacon 
signal at 28.56 GHz has been in continuous operation since March of 1977 
at Wallops Island, Virginia. During the summer of 1977, simultaneovts 
radar and disdrometer measurements were also made and used for predicting 
path attenuation. The best fit values of a and b of the relation k » a Z 
were deduced for each rain period from the raindrop size measurements; 
where k is the attenuation coefficient [dB/km] and Z is the relfectivity 
factor [mm^/m^]. The measured k-Z relations and the simultaneous radar 
reflectivity measurements along the beacon path were Injected into a 
computer program for estimating the path attenuation. Predicted attenu- 
ations when con^ared with the directly measured ones showed generally good 
correlation on a case by case basis and very good agreement statistically. 

The results demonstrate the utility of using radar in conjunction 
with disdrometer measurements for predicting fade events and long term 
fade distributions associated with earth-satellite teleconununi cations. It 

also suggests its use for arriving at sampled fade distributions from which 
prediction criteria associated with variable path angle, frequency, and 
site diversity can be tested. 


vi 



ThI JOHN* HO#KtN( UNIVENItTV 

APPLIED PHYSICS LABORATORY 

Uuxtl. M«RTL/tNO 


PREDICTION OF ATTENUATIC»I OF TOE 28 GHz COMSTAR BEACON SIGNAL 
USING RADAR AND MEASURED RAIN DROP SPECTRA 

1.0 INTRODUCTION 

The ability to use radar as a predictor of rain attenuation has 
been previously tested against actual direct measurements by several 
investigators. McCormick [1] and Strickland [2] correlated their 3 GHz 
radar data with path attenuations at 15.3 GHz; the former investigator 
using a beacon on an aircraft « and the latter using transmissions from the 
ATS-5 satellite. Goldhirsh [3] correlated his 3 GHz radar measurements 
with 13 and 18 GHz measured attenuations using the ATS>6 in an uplink mode. 
Although some statistical aspects were considered by the above investiga- 
tors, none examined how well cumulative fade distributions as derived from 
radar cos^ared with those obtained from direct measuren^nts. More recently, 
Hodge and Austin [4] compared radar measurements at 3 GHz using the McGill 
radar with radiometer-predicted attenuation at 13 GHz and were able to 
arrive at predictive parameters as well as a radar calibration by matching 
radar and radiometer-derived cumulative fade distributions. 

We describe here an experiment performed at Wallops Island, Va., 
using a disdrometer and radar system for estimating the rain attenuation 
of the COMSTAR beacon signal at 28.56 GHz. Drop size distributions 
measured during each rain period enabled the calculation to be made of 
en^irical parameters which were injected into the radar prediction scheme. 

In this way, the dominant drop size distribution for each rain period was 
used. Results using the measured drop size spectra are compared with those 
using the Marshall-Palmer distribution [5]. Individual case histories and 
overall cumulative fade distribution comparisons are presented. The data 
base corresponds to five rain days during the summer of 1977 (June through 
September) and represents 13 attenuation events and 304 minutes of simul- 
taneous radar data. 

2.0 EXPERIMENTAL CONFIGURATION 

2.1 General 

The experimental configuration which is depicted in Fig. 1 con- 
sists of a phase locked loop receiving system operating at 28.56 GHz, an 
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S-band radar (f » 2.84 GHz) located approximately 30 m away» and a system 
of three raingages and two dls drome ters for measuring rain drop size 
distributions located in the immediate vicinity of the receiving antenna. 
Both antennas are fixed and point in the direction of the (X3HSTAR 
geosynchronous satellite [95^W longitude ^ O.l”] with elevation and azimuth 
angles of 41.6* and 210*, respectively. 

In Fig. 2 are shown the geographical location of Wallops Island 
(WI), Va., and the azimith pointing direction. The circled region ^out 
WI represents the radar observation area. 

Tne satellite receiver as well as the disdron^ter-raingage syst«D 
are functioning continuously for arriving at long term statistics (to be 
described in a future report). During selected periods of rai'.i, the radar 
monitors the reflectivity along a parallel path for purposes of testing 
radar methods for predicting the path attenuation at 28.56 GHz. The dis- 
drometer data are used as input in these prediction methods. 

2.2 Receiver System 

The pertinent parameters for the C(»iSTAR receiving system are 
given in Table 1. The system (Fig. 3) consists of an rf front end mounted 
at the antenna location. The IF output of 1.05 GHz feeds a phase lock loop 
receiver located in the radar building. The output voltage rwges from 
0 to 5 VDC and covers the input dynamic range of 32 dB. This output is 
continuously recorded on a calibration strip chart. In addition, the out- 
put signal is coded and recorded on one channel of a four channel tape 
recorder. The tape recorder is automatically actuated whenever the attenu- 
ation exceeds a defined threshold, a raingage bucket tip occurs, or a drop 
impacts on the disdroi^ter sensors. 

Calibration of the receiver system (Fig. 4) is performed by in- 
jecting a stable 1.05 GHz signal at various levels into the front end where 
it is upconverted and routed (switch in up position) through the mixer pre- 
amplifier downconverter. The 1.05 GHz IF output passes through the phase 
locked loop system. Calibrations are performed daily and recorded on tape. 
On several occasions, the satellite beacon signal was used as a source with 
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Figure 1. Experimental configuration at Wallops Island, Virginia. Antennas print in the direction of 
the Satellite at 95" W longitude. 
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Figure 2. Geographic location of test site 
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Table 1 

PERTINENT COMSTAR RECEIVING SYSTEM PARAMETERS 


Antenna Gain 
Beanwidth 
Antenna Diameter 
Free Space Power Received 
ERP 

Path Loss 
Line Loss 
Sky Loss 

Phaselock Hold-In Threshold 
Phaselock Acquisition Threshold 
Dynamic Range (Minimum) 
Predetection 3 dB Bandwidth 
Post Detection 1 dB Bandwidth 
Amplitude Measurement Error 
Noise Figure 


52. 7 dB 
0.4* 

1.83 m (6 ft) 
-106 dBm 

55.8 dBm 
-213.4 dB 
-0.5 dB 
-0.3 dB 
-138 dBm 
-133 dBm 
32 dB 
50-100 Hz 
10-20 Hz 
♦ 1 dB 

18 dB 
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Figure i. Elements of satellite receiving configuration. 









Figure 4, Calibration configuration for the CCMSTAR receiver. 
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the attenuation placed at the rf front end input. These latter calibra- 
tions are within 0.5 dB of the loop calibrations. 

Since the satellite position is maintained fixed through 
stationkeeping to within ^ 0.1®, the receiver antenna is also kept fixed. 
However, an approximate weekly adjustment of the antenna pointing is made 
whenever the satellite assumes a new nodal position. Ihe motion of the 
satellite within its 0.1“ spatial box results in a free space diurnal 
variation of less than 2 dS which is p: dictable to within ^ 0.5 dB. 

The reflection losses due to water on the feed were measured by 
controlled water spray tests and were observed to be on the average of 2 dB 
within an uncertainty of +_ 1 dB using a hydrophobic agent (Fusidox M)*. 

2.3 Radar System 

Ihe operating parameters of the radar (SPANDAR) are listed in 

Table 2. 

During selected periods of rain, the radar reflectivity is 
nwnitored at adjacent range bins along the path where each bin has resolu- 
tion of 150 ra (1 usee pulsewidth). The resulting reflectivity represents 
the average of 128 independent san^les obtained in 0.4 seconds (prf - 320 
Hz) using a frequency diversity transmission scheme where each pulse in the 
prf is shifted by an amount greater than the inverse pulsewidth [6] . The 
integrated data are stored on IBM compatible tape for later analysis on the 
IBM 360/91 computer. 

The electrical parameters of the radar have been directly cali- 
brated [7] and are routinely checked by tracking 6-inch spheres on ballons 
[8] and by monitoring the solar flux from the sum at 2.84 GHz [9]. Rela- 
tive standard gain horn measurements are also periodically performed. The 
accuracy of the absolute calibration of SPANDAR based on on-going measure- 
ments as well as previous component calibrations [7] is well within ^1.5 
dB. 

* 

Silibond Product, Inc., 25 Industrial Way, Wilmington, Massachusetts 
01887 
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Table 2 

OPERATING PARA»CTERS OF SPANDAR 


Radar 

Gain 

Peak Power 

Beamwidth 

Diameter 

Frequency Diversity; 

Center Frequency 
Frequency F.xcursion, Af 
Pulse Length 
Pulse Range Resolution 
prf 

Nunijer of Pulses Integrated 
Averaging Time 

Number of Adjacent Range Gates 
Calibration Error (rms) 


FPS/18 
50.6 dB 
1.0 MW 
0.4“ 

18.3 m (60 ft) 

2.840 GHz 
+ 13.5 MHz 
1 vsec 
150 m 
320 pps 
128 

0.4 sec 
871 

< 2 dB 
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2.4 Disdroaeter-Raingage System 

The raingages used are of the tipping bucket variety and the 
disdrometers which were developed at APL [10,11] are of the iopact type. 

In Fig. 5 is shown a cross section view of the sensor itself. When a drop 
i^acts at tenoinal velocity on the plexiglass aperture, a voltage is 
induced at the piezoelectric transducer and this voltage is prq>ortional to 
the drop momentimi. The transducer output voltage is passed through a pre- 
an;)lifier, shaping circuitry, a biphase encoder, and ultimately to one 
channel of the four channel tape recorder. Calibrations of the disdrometers 
are perfom^d by dropping drops of known diameters through capillaries at 
sufficient heights so that drop terminal velocities are achieved. 


The basic elements of the disdrometer-raingage system are given 
in the block diagram in Fig. 6. Ihe coded outputs of two disdrometers are 
each fed into one channel of a four channel tape recorder. Outputs of 
three raingages are multiplexed with a tin« code (IRIG B Clock) and fed 
into another channel of the tape recorder. Simultaneously the raingage 
and disdrometer rain rates obtained using hard wire circuitry are recorded 
side by side on a chart recording system. A real time comparison of these 
rain rates represents an on-going check on the performance of the disdrometer. 


3.0 ANALYTICAL ASPECTS 

3.1 Radar Near Field Correction 


Since reflectivities are measured at ranges as close as 600 m 
away from SPANDAR, and the far zone distance, r^, is 



(3.1) 


it is necessary to account for near zone effects. In this section we des- 
cribe a method which enabled the measurement of reflectivities to be made 
as close in as 600 m by correcting the range -dependent gain and beamwidth 
in the near field region. 

For a tapered aperture illumination, the gain reduces and the 
beamwidth increases for axial distances smaller than r^. At the shorter 



Figure 5. Cross sectional view of the disdrometer sensor 
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Figure 6. Basic elements of disdrometer-raingage system 
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ranges there results a nonotonic increase in power density (at a rate 
2 2 

smaller than 1/r ) down to 0.1 (2D /X). At closer ranges oscillations 
occur with a diminishing envelope [12]. 

It is apparent from the id}Ove that the radar equation must be 

appropriately modified so as to account for the near field effects. In 

Appendix A we have modified the radar equation using the near field formu- 

2 

lation of Hansen [12] to be applicable down to 0.1 (2D /X) s 600 m. At 
closer ranges the clutter of SPANDAR <k>es not allow any further adjustment 
and we therefore assume the reflectivity to be uniform for the first 600 m 
of range and use the value measured at 600 m. Using the results of 
Appendix A» the near zone power, P(r), is obtained by dij^ding the far 
zone power formulation, P^(r), by the appropriate range dependent correction 
factor, C(r). That is. 


where. 


P(r) * 


P„(r) 

CCr) 





and where. 


r < r 


f 



(3.2) 


(3.3) 


(3.4) 


3(r) 


256 

2 



+ 


128 

2 



(3.5) 


X = 


r 



(3.6) 
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The reflectivity factor in the near zone, Z [dBZ] is related to the 
reflectivity factor using the far zone formulation by the relation, 

Z[dBZ] = Z^[dBZ] + 10 Log^^ C(r) 

The correction factor C(r) may be shown to be the ratio of the far field 
gain to the near field gain. Hiat is. 


CWdB 


10 I^»io 



(3.7) 


(3.8) 


In Fig. 7, C(r)^ is plotted as a function of range from SPANDAR. We may 
deduce from this curve, for example, that at 0.1 (2D^/X) - 600 m from 
SPANDAR, a correction of 4.3 dB must be added to the reflectivity factor 
calculated using the far zone formulation. The correction factor given by 
(3.3) was used in the radar program to arrive at the corrected reflectivity 
factors in the near region as described above. 

3.2 Radar Prediction Formulation 

We here describe the formulations used in arriving at the radar 
and disdroraeter predicted attenuations and subsequent probability levels. 

The radar power recorded in each of the adjacent range bins is 

6 3 

converted into reflectivity factor levels, Z [ran /m ] and these are injected 
in the formulation for path attenuation given by. 


N 

■ 23 

i=l 


k. Ar 
1 


[dB] 


(3.9) 


where 


k. 

i 



[dB/km] 


(3.10) 


NEAR FIELD CORRECTION (dB) 



Figure 7. Near field correction factor as a function of range from SPANDAR. 
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and where 

A|^(t) > the radar predicted attenuation in tise» t [dB] 
r^ > the attenuation inefficient pertaining to i^^ range 
bin [dB/km] 

Z(rj) ■ the reflectivity factor at range bln 

Ar « range resolution Interval (150 a) 

N > nuaber of range bins used in suaaation (see text) 

The Halt, N, in (3.9) is obtained in the following way. Throt^ 
the observation of photographs of RHI's taken at 1/2 hour to 1 hour inter- 
vals, it is noted as to whether the rains are either convective or non- 
convectivej the latter showing a bright band and the foxier showing cellular 

i ** * », 

activity rising to elevations above the zero degree isotherm height. The 
s unnation is taken to a range corresponding to a zero degree isotherm 
height when a bright band appears since the dominant presence of ice above 
this height causes negligible attenuation. On the other hand, the suaaation 
is taken along the entire path when the rain is convective; implying the 
dominant presence of water above the zero degree isotherm height. This 
procedure has been found to give good results both in the previous work [3] 
as well as in the present e;q}eriaent. 

The values of a and b were arrived at by s sibling continuous 30 
second drop size distributions (hereafter referred to as DSD) with the APL 
disdroaeter during each rain period and calculating 

D 

/kaax 

k - y N(D) C^^^(D) dD [dB/km] 

D . 
min 

D 

A fi z 

N(D) D° dD [mm°/m^] 

D . 
min 


(3.11) 


(3.12) 
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where 

k a theoretical attenuation coefficient [dB/km] 

N(D)dlD > nu^er of drops per unit volume between dr<^ sizes 
D and D+dD (DSD) 

3 

■ extinction factor (dB/km) x cm 

6 3 

Z = theoretical reflectivity factor [mm /m j 

D_. , D c the minimum and maximum diameter n^asured, 

min max 

respectively. 

Best fit values of a and b were generated for each rain period; the rain 
period representing a continuous rain event for the day in question. The 
extinction factors for 28.56 GHz, C^j^^(D), were arrived it by inteipolating 
the results of Medhurst [13]. 

4.0 RESULTS 

As mentioned previously, the results presented here correspond to 
the summer of 1977 during which approximately five hours of radar data of 
rain were obtained simultaneous with path attenuation measurements. This 
period enco*^assed five rain days and 13 attenuation events. The disdrometer 
results presented were taken from data acquired over a period which flanked 
the radar measurements. 

In Table 3 are summarized the san^ling days and times for both 

the radar and disdrometer measurements. Also listed are the measured best fit 

parameters a and b (5th and 6th columns) pertaining to the relation (3.10) 

for the indicated rain events. In the 7th column are given values of the 

2 

coefficient of correlation, r , which establishes a measure of the goodness 
of fit; 1 being perfect. In establishing the number of samples and to 
ensure a statistically adequate drop size number, only those spectra were 
selected in which the rain rate exceeded 2.5 irnn/hr. 

4. 1 Empirical Relations from Drop Size Spectra 

Exan^les showing scatter plots of log k versus log Z are given 
in Figs. 8 through 11 for days 157 (June 6), 236 (August 24, afternoon and 
evening rain periods), and 257 (September 14). Each point is representative 
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Table 


SUMIARY OF SANPUNG DAYS AND TIMES A 


Date 


SaiB>ling Time Intervals 

/Jullan\ 

Event 





\ Day / 
1977 

No. 

Radar 

DisdroTOter 

1 

June 6 

B 

21:26-21:55 

21:26-22:45 

(157) 

H 

21:57-22:26 



H 

22:36-22:48 


June 9 

B 

12:44-13:13 

12:44-14:49 

(160) 

B 

13:24-13:53 



B 

14:41-14:50 


Aug 24 


14:56-15:24 

14:20-16:30 

(236) 

8 

15:38-16:07 



9 

16:09-16:25 


Aug 24-25 

10 

23:35-23:59 

23:30-00:52 

(236-237) 

11 

00:05-00:34 


Sept 14 

12 

20:03-20:32 

20:03-20:46 

(257) 

13 

20:33-20:48 


Overall 
Data Base 

304 min 

458 min 

Marshall-Palmer 
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D CORRESPONDING DISDROMETER RESULTS 


Disdroneter Best 

Fit Parameters 

a 

b 


No. of 
Samples 

1.17x10’^ 

0.825 

0.93 

123 

7.36xl0‘^ 

0.^17 

0.75 

180 

1.23x10“^ 

0.492 

0.60 

102 

1.45x10"^ 

0.806 

0.94 

107 

3.47x10’^ 

0.733 

0.95 

53 

1.87x10"^ 

0.775 

0.89 

565 


2.01x10-3 


0.773 
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of a 30-second drop size distribution where the corresponding calculations 
are made using (3.11) and (3.12). The solid line shown in each of the 
figures represent the best fit linear regression line used for that rain 
period. 

In Fig. 12 are shown the best-fit k-Z curves plotted side by side 
for the five rain days. Also shown is the Marshall -Palmer (hereafter 
referred to as M-P) best-fit curve. We note that the distributions for three 
of the five rain periods (June 6 , August 24-25, and September 14) are close 
to the M-P line. It may be noted that when all the drop size distributions 
are considered together, the values of a and b (Table 1) are almost identical 
with those corresponding to the M-P distribution (i.e., a ■ 2.10 x 10~^, 
b = 0.773). Tbis is attributed to the fact the M-P distribution is repre- 
sentative of a large data base for the geographic region of Washington, D. C. 
(mid-Atlantic east coast). On an individual rain period basis the values of 
a and b may be markedly different from the M-P distribution as is indicated 
for the June 9 and August 24 periods. 

Other interesting features of the above results are that the drop 

size distributions for two distinct rain periods occurring for the same day 

may be markedly different as evidenced by the results of August 24. It also 

may be noted from Table 1 that the smaller the coefficient of correlation, 

2 

the greater the deviation from the M-P, and where r exceeds 0.9> the fit 
is close to M-P. 

4.2 Measured and Predicted Attenuation Events 

In Figs. 13 through 17 are given examples of attenuation versus 
tinffi plots of the predicted and measured attenuations corresponding to the 
best-fit regression parameters a and b obtained from the scatter plots of 
Figs. 8 through 11. The ordinate and abscissa in each of these curves 
represent the attenuation in dB and time in GMT, respectively. The solid 
curve represents the measured level, the dashed curve the predicted level 
using (3.9) and the measured empirical parameters, a and b. Also plotted 
is a dot- dashed curve which represents the attenuation obtained using (3.9) 
and a and b corresponding to the M-P distribution. The directly me'.sured 
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Note; for k = aZ 
a = 1.23 X 10-2 
b = 0.492 
r2 = C.60 


Best fit for k = aZ^ 


Log^Q of reflectivity factor, Z [mrr 'm3] 


Figure 9. Scatter plot of log k vs log Z for early afternoon of 
Aug 24, 1977 (day 236) from 14 hr, 21 .dn to ' ' hr, 30 
rain GMT. Solid line is best fit for k » a rela';ion 
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Best fit for k = aZ^ 


Note: for k = aZ 
a = 1.45 X 10-3 
b = 0.806 
r2 = n Od 


Log^Q of reflectivity factor, Z {mm®/m3] 

Figure 10. Scatter plot for log k vs log Z for evening of Aug 24 and 
morning of Aug 25, 1977 (days 236-237) from 23 hrs, 30 min 
to 00 hr, 52 min GMT. Solid line is best fit for k = a 
relation. / 
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Figure 11. Scatter plot of log k vs log Z for Sept 14, 1977 
(day 257 ) from 20 hr, 03 min to 20 hr, 46 min GMT. 
Solid line is best fit for k ® a 2*^ relation. 







" 21:27 21:29 


Figure 13. 


21:31 21:33 21:35 21:37 21:39 21:41 21:43 21:45 21:47 21:49 21:51 21:63 21:56 

Time (GMT) 


Con|>axison of directly measured fade event with DSD-radar predicted and 
M-P-psredicted for June 6, 1977 (day 157) from 21 hr, 26 ndn to 21 hr, 

55 min GMT. 
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Figure 14. Comparison of directly measured fade event with OSD> radar predicted 

and M-P-radar predicted for June 6, 1977 (day 1S7) from 21 hr, 57 min 
to 22 hr, 25 min GMT. 
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Time (GMT) 


Con;>arison of directly measured fade event with DSD-radar 
predicted and M-P-radar predicted for early afternoon of 
Aug 24, 1977 (day 236) from 16 hr, 09 min to 16 hr, 26 min 
C3MT. 


Figure 15. 



Directly measured 
DSD - radar predicted 
MP - radar predicted 


23:36 23:38 23:40 23:42 23:44 23:46 23:48 23:50 23:52 23:54 23:56 23:58 * 24 ; 

Time (GMT) 


Figure X6. Comparison of directly measured fade event with DSD-radar predicted 
^d M-P-radar predicted for evening of Aug 24. 1977 (day 2361 from 
23 hr, 35 min to 23 hr, 59 min GMT. 



20:04 20:06 20:08 20:10 20:12 20:14 20:16 20:18 20:20 20:22 20:24 20:26 20:28 20:30 


Time (GMT) 


Figure 17. Conparison of directly oieasured fade event with DSD- radar predicted 
and M-P-radar predicted for Sept 14, 1977 (day 257) fro« 20 hr, 03 
nin to 20 hr, 32 min GMT. 
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curves are plotted from receiver data saiqtled at 3 second intervals and the 
radar predicted curves are plotted from data sai^les at 10 second intervals. 

The following interesting features may be noted from t!ie above 
comparisons: (1) The receiver dynamic range limits the maximira attenuation 

excursion due to rain to approximately 30 dB beyond which a loss of lock 
condition results. This accounts for the flattened sections of the curves. 

The radar predicted levels do not suffer from this condition and show 
attenuations sometimes in excess of 72 dB. (2) The directly measured and 
predicted curves using DSD and radar measurements show overall good agree- 
^nt although slots of time do exist exhibiting large <teviations. (3) On 
those days in which the measured drop size spectra differed from M-P, the 
radar predicted attenuation using the M-P distribution showed poor ^reement 
with the directly measured attenuation. An exasple of such case is depicted 
in Fig. IS. On the other hand, the DSD predicted results show rather close 
agreement. (4) For June 6 (Figs. 13 and 14), evening of August 24 (Fig. 16), 
and Septend>er 14 (Fig. 17), the predicted attenuations using the measured 
distributions are practically coincident with those of M-P. This, of course, 
should be expected, as the measured values of a and b are similar to those 
of M-P. 

4.3 Measured and Predicted Probability Distributions 

We here compare the cumulative fade distributions of the directly 
measured and radar predicted cases. Such a comparison is given in Fig. 18. 

The vertical scale represents the probability the attenuation exceeds the 
abscissa; the abscissa corresponding to various attenuation levels. The 
solid points, circles, and square points represent, respectively, the 
directly measured, the DSD-radar predicted, and the M-P-radar predicted 
distributions . 

The table in Fig. 18 summarizes quantitatively how well the radar 
predicted fade distributions agree with the measured ones. We note the rms 
and average dB deviations for the radar-DSD case are 1.2 and 1.0 and the 
average ratio of probabilities associated with the predicted and measured 
points is 1.1. For the M-P case, the respective values are 5.6, 5.0, and 1.37. 
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Fade depth (dB) 


Figure 18. Cojiq)arison of cumulative distributions for directly 
nwasured, DSD*radar predicted, and M-P>radar pre- 
dicted cases covering five rain days, 304 minutes 
of simultaneous radar sampling during suaraer of 1977 
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Hie average probability ratio of the two curves was obtained by averaging 
the ratio of the larger probability to the smaller one for eadi dB level vp 
to 28 dB. This gives a measure of the percentage uncertainty of the pre- 
dicted curve relative to the measured one. Hence, on an absolute basis, 
the disdrometer curve shows very good agreement with the directly measured 
whereas the M-P curve shows poorer agreement. Whether this would tend to 
obviate the use of the M-P spectra for such radar predictions is too early 
to say, as an increased data base might bring the M-P radar predicted 
distribution in closer agreement with the directly measured one. 

4.4 Empirical Calibration Adjustment 

An adjustment in the radar calibration was deemed necessary in 
arriving at the reflectivity factors, Z; otherwise, the radar predicted 
attenuation levels consistently underestimate the directly measured values. 
When the reflectivity values were increased by multiplying by an empirically 
determined fixed factor of 3.5 (5.4 dB), both the case-by-case attenuation 
events as well as the cumulative distribution showed good agreement. The 
reason for this required adjustment is not at present understood although 
it might in part be due to a phenomenon described by Rogers [14]. Rogers 
demonstrated that for the case in which the backscatter is measured through 
a logarithmic receiver and there are simultaneously large gradients of reflec- 
tivity in the radar pulse volume, estimated reflectivity levels deviate from 
the true averages by an amount depending upon the gradient value of Z in the 
pulse volume. For example, according to one model, if there is a logarithmic 
variation of Z in the pulse volume with a ratio of maximum to minimum Z of 
6, this could account for a 6 dB underestimation of the reflectivity factor. 

Another possible reason for underestimations of the peak attenu- 
ation values may be due to the con^arative spatial resolutions of the 
receiving and radar systems. The attenuation of the received beacon signal 
takes place essentially through the first Fresnel spheroidal region whose 
diameter is given by 

D(FZ) - 2(r^ - 0.205 r^^^^ [m] 
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where [m] is the axial distarce from the receiving antenna and A [a] is 
the beacon wavelength. The diaaeter of the cross section of the radar beam- 
width (0 ■ 0.4*) is given by, 

Do - 6.98 X 10“^ r^ [a] 

K 0 

and the ratio of these diameters is, 

dW ■ W 

We note that at a distance of 6 km the diameter associated with the radar 
beasMidth is approximately 2.S times the Fresnel zone diaaeter and at 1 km 
they are equivalent. The received beacon signal is therefoz^ expected to 
have greater attenuation peaks than that estimated by the radar because of 
the differences in averaging that takes place across the beams. 

5.0 SUMMARY. DISCUSSION. AND CONCLUSIONS 

Although the overall data base was limited, the DSD and radar 
results showed the following salient features: 

(1) The radar predicted attenuations using the measured rain 
drop spectra showed generally good agreement on a case by 
case basis with the directly measiured levels; although 
sluts of time exist exhibiting deviations. On the other 
hand, extended rain periods exist where the M-P radar pre- 
dicted curves show overall poor agreen^nt. 

(2) The DSD-radar predicted cLumlative fade distributions show 
very good agreement with the directly measured ones whereas 
the M-P results showed poorer agreement. 

2 

(3) When the coefficient of correlation was high (i.e., r > 0.9), 
the measured disdroneter results were in close agreement with 
the M-P spectra; and conversely. 

TVo distinct advantages of using radar results appear in the curves 
shown in Figs. 13, 14, 16, and 17. First, no limitation for examining severe 
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attenuation levels exists in the radar predicted cases as does exist in a 
fixed dynaaic range receiver. Seccmdly, the receiver used could not 
reacquire lock imtil it was ^ to S dB below the leas of lock condition and 
tills rea^uisiti<m ni^t take place only after a frequent search which 
MUld take seveml ainutes. The radar aeasuTMMits do not suffer froa these 
intrinsic receiver characteristics, (hi the other hand, the radar predicted 
aeasureamts do suffer froa other features and these aira: (1) An absolute 

calibration of the radar is necessary, (2) The dre^ size distribution may 
be both spatially and t^orally variable, (3) The attenuation for the 
received satellite signal occurs along a narrow Fresnel ellipsoid whose 
aaxiauB cross section diaaeter is 20 a whereas the radar prediction aeasures 
the average reflectivity Bf the radar beaawidth whicH iriy have ^roxiaately 
3 tines the cross section diaaeter, (4) The radar reflectivities are avezaged 
over a pulse voluae defined by the antenna beaawidth and the range resolutiem 
interval of 150 a within whidi the reflectivities may be hii^ly non-imifora, 
(5) The bean axis of the radar and the receiver systea are not coincident; 
being displaced by 30 a, (6) Near field correction foraulations applied to 
the radar antenna system are based on idealized assuaptions, and (7) The 
pulse voliaes aay cemtain snow or ice near the zero degree isothera height. 

In spite of the nany reasons for non-correlaticm, the radar results 
do correlate remarkable well for the ti>ove example. In fact, the above 
results demonstrate the utility of using radar coupled with disdroneters 
neasurenents for predicting individual fade events as well as long tern fade 
distributions associated with satellite coa«micati<»is through rain. It 
also suggests its use for arriving at saapled fade distributions from which 
prediction criteria associated with varitiile frequencies, path angles, and 
site diversity can be tested [IS]. 

In the absence of DSD spectra, the N-P distribution might represent 
a viable alternative for use with the radar measurements; at least for the 
mid*Atlantic coast geographic xwgiem. An insufficient data base, however, 
exists here to be definitive about this assessmnt. 
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Appendix A 

DERIVATIONS OF THE RADAR EQUATION IN THE NEAR FIELD OF SPANDAR 

The antenna of SPANDAR may be approximated by a circular aperture 

with a tiered illumination having the cross section distribution given by 

where p is the cross section radial distance and D/2 is the 

radius. Using the results of Hansen [12] the gain* G(r)« and beamwidth, 

2 

8(r), in the near region of this aperture (i.e., r ^ 2 D /X) are range 
dependent and are given by. 


0(r) = Y(r) 


G(r) = 


G 

o 

2 

Y 


where 



where 



and where G^ and 9^ are the far field gain and beamwidth, respectively. 

We assume that the radar equation in the near region to have the 
same form as that in the far region with the exception that the beamwidths 
and gains are range dependent. That is, 


P 

r 


( — h ^ 

\l024 Ti Zn 2/V ^ / 


G‘(r) e^(r) 0j^(r) 


(A.l) 
(A. 2) 

(A. 3) 

(A. 4) 


(A.S) 
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where 

c ■ the velocity of light (m/sec) 

= the transmitted power (watts) 

T s the pulse width (sec) 
n ■ the reflectivity (m”^) 

X > the wavelength (m) 
r = the range (n) 

Dividing (A. 5) by the far field formulation, we obtain 



G^(r) e^(r) 6j^(r) 

^ 

G ‘ e 0. 
o eo ho 


where the subscript o denotes the far field values. 

Substituting (A.l) and (A. 2) into (A. 6) 



(r) ^o 



6(r) 


where from (A. 3), 6(r) is given by. 


3(r) 



(A.6) 


(A. 7) 


(A. 8) 


We note from (A. 7) that the far zone power formulation, P^, is 

modified by a correction factor equal to the ratio of the near to the far 

2 

field gains. In Fig. 7, 10 Log^^ of y (r) is plotted down to 600 m from 
SPANDAR. We note that this parameter monotonically increases down to this 
range reaching a value of 4.3 dB. The far zone power formulation must 
therefore be reduced by the indicated levels in the near zone region. 


